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Abstract 

We present a method, based on loop equations, to compute recursively, all the terms 
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1 Introduction 



Random matrix models [33 Ell 1131 EZ1 El E3 have a wide range of applications in 
mathematics and physics where they constitute a major field of activity. They are 
involved in condensed matter physics (quantum chaos localization, crystal 

growths OS ....etc), statistical physics [I31II21IH1E3 ( on a 2d fluctuating surface, also 
called 2d euclidean quantum gravity, linked to conformal field theory), high energy 



physics (string theory [131221, quantum gravity [131123123, QCD [H],...), and they 



are very important in mathematics too: (they seem to be linked to the Riemann con- 
jecture jHSIHIl), they are important in combinatorics, and provide a wide class of 
integrable systems [321E1E3II3- 

In the 80's, random matrix models were introduced as a toy model for zero- 
dimensional string theory and quantum gravity f21EJIH]- 

The free energy of matrix model is conjectured 2 to have a 1/N 2 expansion jlUllTl|ITll2"] 
called topological expansion (N is the size of the matrix): 



That expansion is the main motivation for applications to 2-dimensional quantum 
gravity [HI EI]) because each is the partition function of a statistical physics 
model on a genus h surface. 

The authors of pQ invented an efficient method to compute recursively all the F^'s 
for the 1-matrix model, and they improved it in [2]. 

Their method was generalized in ^H] for the 2-matrix model in the so-called 1-cut 
case. 

Here, we extend the result of [TH], to multicut cases. 

Assuming that the 1/iV 2 expansion exists, the aim of the present work is to give a 
method to compute recursively the terms of the expansion, similar to that of [TH] . 

The 2-matrix model [TIHITT] was first introduced as a model for two-dimensional 
gravity, with matter, and in particular with an Ising field [2D|S]- The diagrammatic 
expansion of the 2-matrix model's partition function is known to generate 2-dimensional 
statistical physics models on a random discrete surface [THIT2"1 E?T] : 



2 There is at the present time no rigorous proof of the existence of the topological expansion; The 
Riemann-Hilbert approach seems to be the best way to prove it as in |17|. The Ricmann-Hilbcrt 
problem for the 2-matrix model has been formulated 6,5,4,30], and seems to be on the verge of being 
solved @|. 
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where the Action is the matter action (like Ising's nearest neighboor spin coupling) plus 
the gravity action (total curvature and cosmological constant) ^3]. The cosmological 
constant couples to the area of the surface, and N (the size of the matrix) couples to 
the total curvature, i.e. the genus of the surface. The large N expansion thus generates 
a genus expansion: 

CO 

F = ^N~ 2h F {h) (1.3) 

h=Q 

where F^ is the partitrion function of the statistical physics model on a random 
surface of fixed genus h. 

p(h) _ y2 \2 e" Action (1.4) 

genus h surfaces matter 

The leading term F^ computed by [3] (along the method invented by and rig- 
orously established by |28j) is the planar contribution. Our goal in this article is to 
compute F^ and present an algorithmic method for computing F^ for h > 1. We 
generalize the method of [TH] . 



1.1 Outline of the article 

• In section El we introduce the definitions and notations, in particular we define 
the 1-loop functions and 2-loop functions, loop-insertion operators, and we write 
the "Master loop equation". 

• In section |3J we observe that, to leading order, the master loop equation is an 
algebraic equation of genus zero, and we study the geometry and the sheet- 
structure of the underlying algebraic curve. 

• In section El we include the previously neglected 1/N 2 term in the loop equation, 
and we compute the 1-loop function Y(x) to next to leading order. Then we 
derive the next to leading order free energy F^ by integrating Y^(x). We also 
discuss how to compute higher order terms. 

• In section |H1 we complete the calculation for the case where the algebraic curve 
has genus one. 

• section |S] is the conclusion. 
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2 The 2-matrix model 



Let N be an integer, V\ and V 2 two polynomials of degrees d\ + 1 and d 2 + 1: 

V 1 (x)=^ +X;|a: fc , V 2 (y) = sfo + £ 1 |y* (2.1) 
fc=i fe=i 

Then let g be an integer (called the genus) choosen between < g < d\d 2 — 1, and let 
e be a g dimensional vector: 

e= (e 1 ,e 2 ,...,e g ) t (2.2) 

We define 

9 

e s +i = 1 -5> ( 2 - 3 ) 
1=1 

All those numbers given, we define the partition function Z and the free energy F 
by the following matrix integral: 

Z = e~ N2p = J dM 1 dM 2 e~ NtI tWHW)-^] (2 4) 

where the integral is over pairs of hermitian matrices Mi and M2 restricted by the 
following condition: 

- the large N limit of the density of eigenvalues of Mi has a support made of g + 1 
disconnected intervals. 

- the filling fraction (integral of the density) in each interval is 6j. 

We assume that the free energy F admits a topological 1/N 2 expansion: 

F = F (0) + _l F (i) + ... + _i_ F W + ... (2 . 5) 

That asumption plays a key role in many areas of physics, in particular quantum 
gravity or string theory fSJE], an d is believe to hold for a wide class of potentials. 
However, the existence of the 1/iV 2 expansion for the 2-matrix model has never been 
proven rigorously (for the one matrix model, it has been established by |17j). 

The goal of this article is to develop a method to compute F^ by recursion on 
h. In particular, we will explicitely compute F^\ Notice that this was already done 
in ^H] in the case g = 0. Notice that F (0) was computed by [HJ. 

Remark: The model can be extended to normal matrices with support on complex 
paths, i.e. the eigenvalues are located along some line in the complex plane, not 
necessarily the real axis. In that case, the potentiasl need not have even degrees and 
positive leading coefficient, the potentials can be arbitrary complex polynomials, and 
the complex paths have to be chosen so that the partition function eq. (|2.4|) makes 
sense. 
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2.1 Definition: resolvents 

We define: 



T M :=l(trM*M<) (2.6) 



The resolvents are formaly 3 defined by: 

Ti rp 
k,0 T T 7" I \ J 0,/ 



in other words: 



fc=0 i=o y 



iV \ x-Mxf K J N\ y-M 2 

We also define: 

:= - W^x) , X(y):=Vl(y)-W 2 (y) (2.9) 

We assume that all the T^ i have a 1/N 2 expansion, and we can write (formaly): 

Y(x) = Y^(x) + + . . . + + . . . (2.10) 

*G/) = * (0) (y) + + • • • + ^ (fc) (y) + • • • (2-ii) 

We will recall below that the leading terms Y^(x) and X^(y) are solutions of 
algebraic equations. Then we will explain how to compute the first subleading term 
Y^(x). We will show that we can compute all Y^> by recursion on h. 

2.2 Other 1-loop functions 

We define the following formal functions: 

°° °° T 1/ 1 1 \ 

rrr ^ V^V^V^ * ^"'t 1 / + 1 ^ Z W \ _ 
U(x, l/):=EEL = n \ tr x-M, ,-M 2 / ^ 

fc=0 j=0 1=0 x v y t i 

d\ d,2 i—1 j—1 

p(x, y ) ■.= EEEE^^;,^ 1 ^'- 1 -^, 

\ j V ; {Ml ) y>(y) - y 2 '(M 2 ) \ 

" iV\ tr x-M 1 y-M 2 J l ^ i4j 

Notice that P(x, y) is a polynomial in x and y of degree d± — 1, d 2 — 1. 
All those functions have a 1/iV 2 expansion. 

3 Formaly means the following: the sums in the RHS arc not necessarily convergent. Wi(x) is 
merely a convenient notation to deal with all T^o at once. 
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2.3 Loop insertion operators 

We define formaly the loop insertion operators: 
d 1 d k d d 



18 y> k d 



dVi(x) ' xdg ^x k+1 dg k ' dV 2 (y) ' ydg Q f^y k+1 dg* k 

These formal definitions actuary mean that for any observable /: 

Of 



(2.15) 



df_ 

dg k 



Res x k ^ T /, : dx 



dVy(x) 

In particular with the free energy, we read from the partition function: 



WAx) 







dVt(x) 

2.4 2-loop functions 

We define the following functions: 



W 2 {y) 



d 



dv 2 {y) 



(2.16) 



(2.17) 



Q(x; x 1 ) :- 
n(y;x') :-- 







dVi{x' 




-WAx) = - ( tr 



1 



tr 



1 



x - Mi x' - Mi 



8V\ far 



^W 2 {y) 



1 1 
tr — — tr 



U(x,y;x') 







dV x {x 



7\U(x, y) 



x' -M x y- M 2 



tr 1 YM ~ ^(M 2 ) fa I 



x - Mi y-M 2 



x' - Mi 



:2.i8i 



;2.i9) 



(2.20) 



2.5 Master loop equation 

It is shown in |18|ll9ll2"U] that we have the following system of equations, called the 
"master loop equations" 



(2.21) 
(2.22) 

(2.23) 



E(x,Y(x)) = —U(x,Y(x);x) 



U(x,y) = x- V 2 \y) + 



U(x,y;x') 



E(x,y) 1 U(x,y;x) 



y - Y(x) N 2 y- Y(x) 
d 



dVi (x 



where E(x,y) is a polynomial in x (degree d\ + 1) and y (degree d 2 + 1): 



E(x, y) := (V((x) - y) (V 2 \y) - x) - P(x, y) + 1 



(2.24) 



and where P(x,y) was defined in eq. (|2.14|) . 

We will solve that system of equations, order by order in 1/N 2 : 




(2.25) 



Then, the large N expansion of the free energy is obtained from eq. (j2.17|) : 







F {h) = 5 h , V{(x)-Y^(x) 



(2.26) 



dVi(x) 



3 leading order: algebraic geometry 



To leading order in 1/N 2 , the master loop equation is an algebraic equation for Y^°\x): 



E(°\x,Y(°\x)) = 



(3.1) 



where 



EW(x,y) = (V{(x) - y){V±(y) - x) - P®\x,y) + 1 



(3.2) 



and P<®(x,y) is a polynomial of degree {d\ — 1, d 2 — 1) whose coefficient of x dl 1 y d2 1 
is known and is equal to g^+ig^+i (from eq. ()2.14|) ). 

Before going to order 1/N 2 and higher, we need to introduce some concepts of 
algebraic geometry, and study the geometry of the above algebraic equation. 

3.1 Determination of PW(x,y) 

So far, we don't know the other d\d% — 1 unknown coefficients of P^°'(x,y). They are 
determined by the following requirements: 

- E(°\x,y) = is a genus g algebraic curve. If g is less than the maximal genus 4 
we get d\d% — 1 — g constraints on the coefficients of p(°\x, y). 

- The contour integrals of Y ( -°\x)dx along A-cycles are: 



where Ai,B iy i = l,...,g is a canonical basis of irreducible cycles on the algebraic 
curve. We get g independent equations for the coefficients of P^(x,y). 

Therefore we can determine the polynomial P^ (x, y) (and thus the polynomial 



4 Thc maximal genus is d\d2 — 1. It can be computed by various methods. In particular, see |32) 
for the Newton's polygons method 




(3.3) 



(x, y)) completely. 
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A physical picture is that the support of the large N average density of eigenvalues 
of the matrix M\ is made of g + 1 intervals [a*, 6j], i — 1, . . . , g + 1, and for each i, A% 
is a contour which encloses [a», 6,] in the trigonometric direction, and which does not 
enclose the other aj or 6j with j ^ i. eq. (|3.3j) means that the interval [a«, bj\ contains a 
proportion of the total number of eigenvalues, this why is called a filling fraction. 

In practice, we don't have a closed expression for P (0) (x,y) function of the 
coefficients of V% and V2 and the filling fractions £j. The converse is easier: given a 
genus g algebraic curve, we can determine V%, V 2 and the e^S. 

3.2 Algebraic geometry 

Before proceeding, we need to study the geometry of our algebraic curve. 

Let us call £ the algebraic curve, and we consider that an abstract point p G £ is 
a pair of complex numbers p = (x,y) such that E^(x,y) = 0. Thus, x = X{p) and 
y = y{p) are complex- valued functions on the curve. 

Here, we sumarize some well known properties of Riemann surfaces and theta- 
functions. We refer the interested reader to textbooks [23121] for proofs and comple- 
ments. 

3.3 Sheet structure 

The function Y^ \x) (resp. X^°\y)) is multivalued, it takes d 2 + l (resp. di + 1) values, 
which we note: 

(F (x),ri(x),...,F d2 (a;)) (resp. (X (y), X 1 (y), . . . , X dl (y)) ) (3.4) 
The zero th sheet is called the physical sheet, it is the one such that (from eq. (J2.9|) ): 

Y (x) ~ V{{x) - - + 0(l/x 2 ) (resp. X (y) ~ V 2 \y) -- + 0(l/y 2 )) (3.5) 

x^oo x y^co y 

An x-sheet (resp. y-sheet) is a domain of £ on which the function X (resp. y) 
is a bijection with C U {+00}. The curve £ is thus decomposed into d 2 + 1 x-sheets 
(resp. d\ + 1 y-sheets). The decomposition is not unique, and a cannonical possible 
decomposition will be given below. 

For each x (resp. y), there are exactly d 2 + 1 (resp. d\ + 1) points on £, one in each 
sheet, such that: 

X[p) = x ^p = Pk{x) k = 0,...,d 2 
(resp. y(p) = y <-> p = p k (y) k = 0,...,d 1 ) (3.6) 

And thus: 

Y k (x) = y( Pk (x)) (resp. X k (y) = X(p k (y)) ) (3.7) 
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3.4 Points at oo, poles of X and y 

In particular, in each x-sheet, there is a point p such that X(p) = oo. It can be seen 
from eq. (|3.5|) . that there are exactly two such points on £. We define p± such that: 



p + (resp. pJ) is in the x-physical sheet (resp. y-physical sheet), while p_ (resp. p + ) is 
at the intersection of the other d 2 x-sheets (resp. d\ y-sheets). 

3.5 Endpoints and cuts 

The x-endpoints (resp. y-endpoints) correspond to singularities of Y(x) (resp. X(y)), 
i.e. they are such that dY/dx = oo (resp. dX/dy = oo), i.e. they are the zeroes of 
dX{p) (resp. dy{p)). There are d 2 + I + 2g (resp. d\ + 1 + 2g) such endpoints: 

dX(p) = <-> p G {ei, e 2 , • • • , e d2+ i +2g } , dy(p) = <-> p£ {ei, . . . , e dl+ i +2g } (3.9) 

the endpoints are such that 3A; 7^ /, Pk(x) = pi(x) (resp. pk{y) = Pi{y)), i-e. they are 
at the intersection of two sheets. 

3.5.1 Critical points 

In a generic situation, all the endpoints are distinct. If V\, V2 and are chosen so that 
some endpoints coincide, we say that we are at a critical point. Imagine that e is such 
a multiple endpoint, near which dX has a zero of degree q — 1 and dy has a zero of 
degree p — 1, then: 



which is a typical critical behaviour of a (p, g) conformal minimal model. 

From now on, we assume that we are in a generic situation, i.e. all the endpoints 
are distinct. 

3.5.2 Cuts 

The cuts are the contours which border the sheets. Like the sheets, they are not 
uniquely defined, there is some arbitrariness. 

A canonical choice for the cuts is the following: the cuts are the sets of p G £ such 
that 




and 




(3.8) 



Y(x) - y{e) ~ O ((x - X{e)f lq ) 



(3.10) 




v 



(3.11) 
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3.6 Irreducible cycles 

We have already introduced a basis of irreducible cycles, Ai, Bi,i — 1, . . . , g, such that: 

Ap|^ < S >.i ( 3 - 12 ) 

Moreover, we assume that the A-cycles are cuts, and that the A and B cycles do not 
intersect a line L which joins p + and p_. 
We have: 

^- I y{p)dX(p) = €i (3.13) 

and we define: 

Tj := j> y(p)dX(p) (3.14) 
We then define the period-matrix r by: 



1 OVi , 

T » ■= ^ < X15 > 

Remark: 

It is proven in the appendix |H] that: 

dF (o) ! d 2 F (o) 

dej 2m de^j 

where is computed from eq. (|2.17j) . Notice that r is symmetric. 

3.7 Holomorphic differentials 

We define the following differential one-forms: 

dUi(p):=^--^(y( P )dX(p)) (3.17) 

they are holomorphic. Indeed, the pole of y(p)dX{p) = dVi(X(p)) — W\{X (p))dX (p) 
at p = p + is independent of (because V\ is independent of 6j and KesWi(x)dx = 1 is 
independent of therefore dui{p) has no pole at p+. By the same argument, dui(p) 
has no pole at p_, and dv,i(p) is holomorphic. Moreover we have (from eq. ()3.13|) and 
eq. (EH): 

duj(p) = 5i.j , ® duj(p) = T Lj (3.18) 

Remark: 

Anticipating a little bit, it follows from eq. (|3.14jl . eq. ()4.27|) and eq. 1)4.30)) . that: 

1 <9r 

= -^mWW) mp) (3 ' 19) 
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3.8 Abelian differential of the third kind 



On the Riemann surface £, there exists a unique abelian differential of the third kind 
dS, with two simple poles at p — p±, such that: 

Res dS = -1 = -Res dS and Vz / dS = (3.20) 



p+ v- 



We choose an arbitrary point p Q G £ , which does not belong to any cut or any 
irreducible cycle, and we choose a line L joining p + to p-, which does not intersect any 
cycle and does not contain p . Then we define the following functions on £\(UiAi Uj 
Bt UL): 

S(p):= f dS , A(p) :=expS(p) (3.21) 

where the line of integration does not intersect any cycle neither the line L (notice 
that the integral around L vanishes because it encloses p + and p_ which have opposite 
residues). We have: 

S(p) has logarithmic singularities near p + and p_, and is discontinuous along L, 
the discontinuity is: 

5S(p) = 2in peL (3.23) 

S(p) is continuous along the S-cycles, and discontinuous along the A-cycles, the 
discontinuity is: 

5S(p) =Vi P^A , Vi '■= <p dS = Ui(p+) - m(p-) (3.24) 

A(p) has no discontinuity along L, it has a simple pole at p + , and a simple zero at 
p_, therefore the following quantities are well defined: 

7 := lim X(p)/A(p) , 7 := lim y(p)A(p) (3.25) 

Remark: By an appropriate choice of po, it should be possible to have 7 = 7, 
however, we will not make that asumption. 

4 2-loop functions and the Bargmann kernel 

Consider the 2-loop function defined in eq. (|2.18j) : 

rw dWAx) 1 dY(x) I 1 1 \ , AnnS 

Q(x: x r ) = - x = 7 ttt-t^ — 7-4 = -(tr tr > (4.26) 

1 ' ; (x f ) (x-x>) 2 dV x {x') \ x-Mt ^-Mi/ conn 1 ; 
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and define the bilinear differential (where x — X (p) and x' = X(p')): 



B(p,p') : 



dY(x) 
dVjx>) 



doc doc 



(4.27) 



Q(x; x') has the following properties: 

• Q(x; x') = Q(x'; x) is symmetric. 

• since Y(x) has square root singularities near the endpoints e^, Q(x; x') has in- 
verse square root singularities near the endpoints (i.e. simple poles in p — > e k ). 
Therefore B(p,p') has no pole in p = e^. 

• since Wi(x) behaves like - + 0(l/x 2 ) in the physical sheet, i.e. when x — > p + , 
we must have: Q(x; x') ~ 0(l/x 2 ). In particular B(p,p') is finite when p — > p + . 

• since F(x) behaves like V 2 '(y(x)) — + 0(1/F 2 (x)) ~ x when x — > p_, we 
must have that B(p,p') is finite when p — >■ p_. 

• f2(x; x') has no pole at x = x' in the same sheet (i.e. when p = p'). This implies 
that B(p,p') ~ (X (p) — ^(p'^^dx^cfa^p') when p — > p'. 

• since Y(x) satisfies eq. ()3.3|) and e is independent on V\, we must have: 



This allows to determine Q(x; x'). Indeed, B(p,p') is a meromorphic bilinear differ- 
ential on £ , with only one normalized double pole at p = p', and normalized A-cycles, 
therefore B(p,p') is the Bargmann kernel, i.e. the unique meromorphic bilinear differ- 
ential with such properties. 

It can be written (see appendix IHJ) : 




(4.28) 



B(p,p') = didj ln#(u(p) — u{p') — z) dui(p)duj(p') 



(4.29) 



It has the property that: 




(4.30) 



Notice that in eq. (|4.27J) . the derivative is taken at fixed x = X(p). 



11 



5 1/N 2 Expansion 



We are now interested in the 1/N 2 expansion of the free energy and loop functions: 

F = + -LfW + . . . , Y(x) = YW(x) + -^Y^(x) + ... (5.31) 
iV z iV z 

where 

^( ft )=«(x)-rW(,) (5.32) 

So far, we have explained how to compute Y ( -°\x). Once Y^°\x) is known, can in 
principle be computed from eq. ()5.32|) . this has been done in [3]. 

Our goal is to compute Y^(x), FW, and then define a recursive procedure to 
compute Y^ h \x) and F^ for all h > 1. 

5.1 1/iV 2 term 

First we expand the polynomial E(x,y): 

E(x, y) = (x, y) + (x,y) + ... (5.33) 

where E^(x, y) = -PW(z, y) is a polynomial of degree (c?i — 1, <i 2 — 1) whose coefficient 
of x dl ~ 1 y d2 ~ 1 vanishes. And we write similar expansions for all other loop functions, in 
particular U(x,y) and U(x,y;x'). 

U(x, y) = U<® (x, y) + ±17® (x,y) + ... (5.34) 

U(x, y;x') = U<® (x, y;x') + -^[/« (x, y; x') + . . . (5.35) 
Then we expand eq. ()2.21|) to order 1/N 2 : 

E (1) (x,F (0) (x)) + F (1) (x)4 0) (x,F {0) (x)) = ?7 (0) (x,y (0) (x);x) (5.36) 

y(1) = fW(x,yf°)(i)) + [/("'(x,yf»>(x);x) 

4 0) (x,y(°)(x)) 

So far, the polynomial pW(j, y) is unknown, i.e. we have d\d,2 — 1 unknown coefficients. 

We expect that order by order in the 1/N 2 expansion, the resolvent W\{x) = 
V{(x) — Y(x) has no singularities appart from the endpoints, so we require that Y^ 
has singularities only at the endpoints. 

The condition that Y^ has singularities only at the endpoints, implies that in 
eq. (|5.37|) . the poles at the zeroes of Ey°\x, Y^(x)) which are not endpoints should 
cancel. Since there are d\d,2 — 1 such points, we can determine PW(x, y), and thus we 
can determine Y^'(x). In other words, is determined by the condition that Y^ 1 ' 
has singularities only at the endpoints. 
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5.2 The function U(x,y;x') to leading order 

Consider x in the physical sheet, so that Y^(x) = Yq(x). From eq. (|2.22j) we have: 

9E(°>(i,y) 



Notice that 



dUM(x,y) 
dV^x') 



W dY^(x) E^(x,y) 



y-Y(°)(x) dVx{x') {y-Y(°){x)y 



E(°\x,y) = -g* d2+1 l[(y-Y k (x)) 



(5.38) 



(5.39) 



fc=0 



(indeed, both sides are polynomials in y with the same degree, the same zeroes and 
the same leading term). Therefore: 



dE^(x,y) 
dV^x') 



d 2 



-E^( Xl y)J2 



dYJx) 



1 



and thus: 



£7(0) ( 



k=0 



/\ E(°\x,y)^dY k 



dViix') y-Y k (x) 



[x 



y-Y (x) t^dV^x') y-Y k (x) 

Notice that we have considerably simplified the derivation given in 
In particular, when y = Y (x) and x' = x, we have: 



(5.40) 



(5.41) 



U^(x,Y (x);x) = E^(x,Y (x)) £ 



dYJx) 



^dV x {x) Y (x)-Y k (x) 



(5.42) 



5.3 yd) 

Using eq. ()5.37|) and eq. ()5.42|) we have: 



yd), 



x) 



d 2 

E 



dY k (x) 



p(V( Xj Y(°\x)) 

4 0) (x,yW(x)) ' ^m(x) Y (x)-Y k (x) 



(5.43) 



and Y^ 1 ' has poles (of degree up to 5) only at the endpoints. In other words, Y^'(x)dx 
is a one form, with poles only at the endpoints (no pole near p + and pJ). 

5.3.1 Behaviour near the endpoints 

Recall that the endpoints are the zeroes of dX(e) = 0. If p is near an endpoint e k , 
there exists a unique (because we have assumed that the potentials are generic) p' such 
that X{p') = X(p) and p' is near e k . 
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We have 



YW(x(p))dX(p) = y{p) [ W) + O(l) whence, (5.44) 

where B(p,p') is the Bargmann kernel. 
This can also be written: 

rW(*M)^M = Re. ( ^ _ ^ff*^ _ yw) + Q (1) whence* (5.45) 
By adding only 0(1) quantities, we arrive at: 

^ 35 Wp) - *%&) - wi + 0(1) when ^ ei (5 ' 46) 

Since that quantity is symmetric in x and y, we have: 



Y^{x{p))dx{p) + x { - l \y{p))dy{ P ) = Res 5(p ' p,) 



+ Cidui(p) 

i=l 



(5.47) 

where Cj are some constants. Indeed, the difference between the LHS and RHS has no 
pole, it is a holomorphic one-form. 

5.3.2 local coordinate near an endpoint 

Consider that z{p) is a local coordinate near an endpoint e k , we have: 



X(p) = X(e k ) + yA"'(e fe ) + jX>"(e k ) + ^X IV (e k ) + ... (5.48) 



y(p) = y(e k ) + zy'(e k ) + ^y"(e k ) + jy"'(e k ) + ... (5.49) 



B(p,p f ) = (j^f + + • • •) dzdz ' 



(5.50) 



where S(p) is the projective connection. 
X(p') = X(p) implies: 

z' = -z{\ + r k z + r 2 k z 2 + (2r 3 k + t k )z 3 + . . .) (5.51) 

where 

1 X"'(e k ) l X IV (e k ) 1 X v (e k ) 

n ~ 3 X"(e k ) ' Sh ~ 6 *"(e fc ) ' * ~ 60 *"(e fc ) TkSk 1 j 
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That gives: 



dX(p) _ 1 _ 4 X"'(e k ) 



-z 



-3 



cfe 8Ar w (e fc )y( ejfe ) 24A"' 2 (e fc )y(e fc ) 

*'"(e fc ) 2 _ * IV (e fc ) A""(e fc ) y"(e fc ) _ Y"(e k ) 
X"(e k )i X"(e k ) ~T~ X"(e k ) y(e k ) y(e k ) 2 

As 



S(e k ) , 



48^"(e fc )y(e^ 
z~ 2 + 0(l) 



12A"'(e fc )y(e fc 
(5.53) 

This is in principle sufficient to determine K^. 

5.4 Free energy 

Then, we want to find the free energy such that: 

F«(x) = -^— ^ (5.54) 
aVi (x) 

In this purpose, we have to compute the derivatives of various quantities with respect 
to Vi(x), and in particular, how the theta-function parametrization changes with the 
potential V\. 

We conjecture: 

F(1) = -^ ln II F '( e *) ( 5 - 55 ) 

i 

where 

5.5 Higher orders 

Imagine we already know all quantities up to order ft, — 1, and write eq. (j2.21|) to order 
ft: 

h ( 3 

j=0 \ k=0 

= ^2 N -2h+2j U {h- J ) f J2N- 2k Y^{x);x J +0{N~ 2h ) (5.57) 

3=0 \ k=0 J 

The only unknown quantities in that equation are: F^'(s) and E^ h \x, Y^(x)). The 
polynomial E^ h \x,y) must be chosen such that Y^(x) has no other singularities than 
the endpoints, and is completely determined by this consition. That allows to find 
as well as E^ h \x,y) and U^ h '(x,y) to order ft. 
The procedure can be repeated recursively to find Y^ h '(x) to any order. 
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6 Genus 1 case 



Let us recall that the case g = was done in ^H]- The case g — 1 is treated in this 
section. 

We require that E(x, y) = be a genus one algebraic curve. Therefore, there must 
exist an elliptic uniformization. We choose it of the following form: 



x = X{u) = 7 



6(u - 1100)6(11 + Uoo ) d * Hi °( u oo - o-i(O)) 



■ ' '0(t I -M oo )*«(!< + « oo )n,9(l'co + <>i(0)) k ; 

and we denote r the modulus. Here, 6 denotes 9±, i.e. the prime form for genus 1. A 
definition of the ^-function and its properties can be found in appendix |HJ 
We must have: 

£Vi(0) = (4-lKc , ^2 o-i(O) = (di - 1)moo (6.2) 

i i 

All this means that for every (x, y) which satisfy E(x, y) = 0, there exists at least 
one u (in the fundamental paralellogram of sides 1, r) such that x = X{u) and y = y{u). 
An alternative parametrization is: 

X{u) = 7 ^y^(u - Uoo ) - Z(u + Uoo )) + A + J2 ^V ( ^ 2) ( M + ««,) (6.3) 



k=2 



di 



y{u) = -7-L^(Z(« + Uoo ) - Z{u - Uoo )) + Ao + JTZ^^iu - ««,) (6-4) 

^ > k=2 

where Z is the Zeta-function, i.e. the log-derivative of 6\, and is the Weierstrass 
function, i.e. = —Z'. 



We note the inverse functions: 

x = X(s) «-> s = a(x) , ?/ = y(s) «-> s = a(y) (6.5) 

The functions cr(x) and er(t/) are multivalued, we will discuss their sheet structure 
below. The functions Y(x) and A(y) are: 

Y(x) = y(a(x)) , = (6.6) 

They are multivalued too, and their sheet structure will be discussed below. 
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6.1 The parameters 

Our parametrization depends on d\+d 2 +6 parameters which are: 0fc(O) {k = 1, . . . , d 2 ), 
^j(O) (j — 1) ■ • • > 7; 7> M oo and r. The condition eq. (|6.2j) means that only di+d 2 +4 
of them are independent. 
Equations eq. ()3.5|) read: 



f DC' 



<//, = T^ ^ ^7^ fc = l,...,d a + l (6i 



and 

1 = J- <f dsy(s)X'(s) = ^-f dsX(s)y'(s) (6.9) 

Where the contour of integrations are small cylces around ±Woo. 
And eq. (|3.3|) reads: 

2me = [ ds y(s)X'(s) (6.10) 
Jo 

We thus have di + d 2 + 4 equations, therefore we can, in principle, determine all the 
parameters. 

In principle, it should be possible to revert these formula, and compute the <t(0)'s 
and cr(0)'s as functions of the coupling constants g and g. This can be done at least 
numerically. 

Remark: Another point of view is interesting too: once eq. ()3.3|) and eq. ()3.5|) are taken 
into account, we have di + d 2 + 2 independent parameters. This is precisely the number 
of coefficients of the potentials V\ and V 2 . We can consider that the parameters are 
merely a reparametrization of the coefficients of the potentials, according to 
eq. flnZj) and eq. (l6~£l) . 

6.2 endpoints and cuts 

The endpoints in the x-plane (resp. y-plane), i.e. the singularities of Y(x) (resp. X(y)) 
are such that: 

X'(s)=0 (resp. y (s) = 0) (6.11) 
There are d 2 + 3 (resp. d\ + 3) such endpoints: 

X'(s) = «-> s E {e 1 ,e 2 ,. . . ,e d2+3 } , y'(s) = «-> s e {e u . . . , e dl+3 } 

(6.12) 
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We can write: 

{) ~ T rii^oc-e i ) 6(s - Uoo ye(s + Uoo y^ { - ) 

7 n^(^oo + e,) ^(s-n 00 )*+^( S + Moo ) 2 l ' ) 

Since A"(s) an d y'( s ) are elliptical functions, we must have: 

J^e* = -(d 2 - 1)^ , J^gj = (di - l)«oo (6.15) 

i i 

and since Af' (s) and y( s ) are the derivatives of elliptical functions, we must have: 
Y,Z{u O0 -e i ) = {d 2 + l)Z{2u oo ) , ^ Z( Moo + et) = (d 1 + l)Z(2u 0O ) (6.16) 

i i 

Notice that eq. ()6.7|) for k = d\ + 1 and eq. ()6.8|) for k = d 2 + 1 imply: 

7 = -di9d 1+ il YT~oi =~\ ' T = - "2</da+i7 77-77 ; V l b - 17 ) 

lli^(«oo - ei) lI i 6'(«oo + e,) 

6.3 2-loop functions and the Bargmann kernel 

Consider the 2-loop function defined in eq. (|2.18j) : 

_ dWxjx) _ 1 dY(x) I 1 1 \ 

(X; X } " dV^) ~ {x-x'Y ~ dVM) ~~\ ^M~i ^ I conn ( } 

f2(x; x') has the following properties: 

• Q(x; x 1 ) = Q(x'; x) is symmetric. 

• since Y(x) has square root singularities near the endpoints Af(efc), f2(x; x') has 
inverse square root singularities near the endpoints (i.e. simple poles in cr(x) — e^). 
Therefore A''(s)f2(A'(s); x') is finite when s — > e^. 

• since Wi(x) behaves like - + 0(1/ x 2 ) when x — > 00 (i.e. a(x) — > +M00)) we 
must have: f2(x;x') ~ 0(1/ x 2 ). In particular Af'(s)0(A'(s); x') is finite when 

S — > +Moo. 

• since F(x) behaves like V^(K(x)) — -p^y + 0(1/F 2 (x)) ~ x when cr(x) — > —Uoo, 
we must have that X'(s)Q(X(s); x') is finite when s — > — Woo- 



f2(x; x') has no pole at x = x' when x and x' are in the same sheet, i.e. when 

dY(x) 1 
dVi(x') ~ (x-x') 2 

1? 



cr(x) = cr(x'). This implies that ~ < x - x >) 2 wnen ° r ( x ) a ( x ')- 



• since Y(x) satisfies eq. ()3.3j) and e is independent on V\, we must have: 

/ n(X(s);x')ds = (6.19) 
Jo 

This allows to determine Q(x;x'). Indeed, the function X , (s)X'(u) qvi(x\u)) * s an 
elliptical function of s, with only one double pole at s = u. Therefore (see appendix 
EJ): 

x '^ x '^Wm = ^- u)+c (6 - 20) 

where <fi is the Weierstrass function, and where the constant C must be equal to zero 
in order to satisfy eq. (|4.28|) . 

We recognize the Bargmann kernel: 



(6.21) 




In a similar fashion, we find that the function: 



= ^4 = - (* tr } (6.22) 

vy ' ; dV^x') \ y-M 2 x'-M 1 / 



is: 




(6.23) 

We are now equipped to compute the next to leading order functions... 

6.4 Computation of 

We have (see eq. (fQ3|) ): 

Y (i) (r) _ P (1) (x,Y (x)) 0(cr fc (x) -a (x)) 

lJ ^(x,F (x)) ^A"K(a;))^( ( 7o(a;))(yo(x)-n(x)) l ' J 

and we require that KW(s)'s only poles (of degree up to 5) are the endpoints Af(efc). 
Note that X' (s)Y^ (X (s)) has no pole at s = ±Woo> therefore we may write: 



, d2+3 

y( 1 )( < V( a )) = __ Vi/( S -e,)+B^(s-eO+^(s-^)+^(«-^ (6.25) 
X'(s< 



k=l 
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The coefficients A k , B k , C k , D k are determined by matching the poles in eq. (|6.24j) 

(P^ doesnot contribute to them because E y (x,Y (x)) has only single poles at the 
endpoints). Below, we will find that D k = 0. 
Let k G [1, d 2 + 3], and choose s close to e k : 

s = e k + e (6.26) 

there must exist s (unique because the potentials are non-critical) such that X(s) = 
X(s) and s is close to e k (s is the cr k (x) in eq. (|6.24j) ): 

S = e k - V , v = 0(e) (6.27) 

By solving X(s) = X{s) order by order in e we get: 

r] = \e , X = l+r k e + r 2 k e 2 + (2rl + t k )e 3 + 0(e 4 ) (6.28) 

where 

1 X"'{e k ) 1 X IV (e k ) . 1 X v (e k ) 

T k = - , Si. = ; — , tu = : : T k S k 6.29 

3 X"(e k ) 6 X"(e k ) 60 X"(e k ) V ' 
From eq. ()6.24j) we must have: 

eAt - 2Bte + ^ + = - x , (et _j { t { :::?_ nek _ v)) + ^ 

We note the 3rd degree polynomial: 

P k (e) = 6A k -2B k e + C k e 2 + D k e 3 (6.31) 



i.e. 



P k (e) = ((l + A)- 2 -Cie 2 2 ) 

(-XX"(e k ) + e^-X"\e k ) - e 2 ^-X IV (e k ) + e ^X v {e k )y l 
2 6 24 

((1 + \)y\e k ) + 5(1 - X 2 )y"(e k ) + f (1 + A 3 )y"(e fc ) 
z 6 

+^(l-AV(e fc ))- 1 + 0(6 4 ) 

= -(^"(e^yCefc))" 1 (l + A)-^- 1 
3 5 
(1 - -eXr k + e 2 X 2 s k - -e 3 (t fc + r^))" 1 

n + e n A) r(efe) + e2 n a + a^^V 1 

(1+ 2 (1 - A) y^y + "6- (1 - A+A) 3^y } 
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It is easy to see that t k as well as y IV (e k ) disappear, and that D k = 0. 
After a straightforward calculation, one finds: 



- 8X"(e k )y'(e k )(QA k - 2B k e + C k e 2 



l 1 X m \e k ) 
S e X"(e k ) 



1 „ fX"\e k f X IV (e k ) 



6^ \X"(e k f X"(e k ) 
X"'{e k )y"{e k ) y"'{e k ) 



(6.33) 



X"{e k )y>{e k ) y>(e k ) 



24Ci 



After substitution into eq. (|6.25j) . we find the genus one correction to the resolvent: 



4SX'{s)X"{e k )y(e k ) 
, ST d 2+ 3 X"'{e k ) 



[s - e k 

it 



s - e k 



~r Z^fc=l 

Ed2+3 
k=l 



48X>(8)X"2(e k )y'(e k y 
X "'(, k f x IV (e k ) x>"(e k ) y" (efc) y"'(e k ) 

d 2 +3 A"'(e fc )2 X"{e k ) ^IFH^J y> (e k ) y'(e fc ) 



mX'{s)X"{e k )y>(e k ) 



[s - e k 



2X'(s)X"(e k )y>(e k y 



(6.34) 



That function represents the partition function of a statistical physics model on a genus 
one surface with one + boundary. 
Notice that 



de 



rd) = E 



d2+3 X"\e k f _ X' v (e k ) X"'(e k ) y"(e fc ) _ y"'(e k ) 
X"{e k Y X"(e k ) "T #»(e fc ) y(e fc ) y(e fc ) 



24Ci 



fc=i 



48A"'( ejfe )y(e fc ) 



(6.35) 



6.5 The free energy 

We are now going to find the free energy such that: 

dF (i) 



r (1) (x) 



dVAx) 



(6.36) 



conjecture: 



1 / d 2 +3di+3 / 



0(ej - e^OQu^) 
9(ci - u 00 )6'(e j + Uoo) 

question: how to extend that conjecture to genus g > 1 ??? 



(6.37) 
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7 Variations with respect to the potentials 

7.0.1 Variations with u fixed 



Now, we would like to compute: 

dX(u) 



X{u) := X'(s) 



and y(u) := X\s) 



dy{u) 



dv l {x{ s )) ' " v y - v ' dV^Xis)) 

with u fixed (we will not write the s dependence for lisibility in this section). 



(7.38) 



We have: 
dy(u) 



dY(X(u)) 



dV^Xis)) 
which implies: 



dV r {X{s)) 



dY(x) 



dV^Xis)) 



dX(u) 



(7.39) 



x{u)y'{u) - y{u)x'{u) = 4>{s - u) 



In particular at u = e« we have: 

X(e t ) 

We set: 

ati : 

Moreover, we have: 

X(u+1) _ X{u) 
X'(u + 1) ~ X'{u) 



s - d) 



0(g - gj) 
X"(e t )y'( ei ) 

X(u + t) _ X{u) 

X'(u + t) ~ X'(u) 



(7.40) 

(7.41) 
(7.42) 

(7.43) 



where f := X'(s) dVl ^( s )) ■ That means that the function: 

<2 2 +3 



X(u) \ 



has no poles and satisfies: 

/(«+!)=/(«) 



/(« + r) = /(«)-f + 2i7r5] 



(7.44) 



(7.45) 



/'(u) is an elliptical function with no pole, therefore it is a constant, and f(u) is a 
constant too. 

That allows to write: 



X(u) = X\u) 



cfc+3 

C + otiZ(u — 
i=i 



(7.46) 
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that implies: 



d 2 +3 

2iir a. 



(7.47) 



and C is a constant which is determined by the behaviours near ±1*00- Indeed, consider 
X'(u)/X'(u), you get: 



2u' n 



and 



(d 2 + l)uc 



u — u r 



d 2 + 1 



U — Ur 



C + y^ j a i Z{-u 00 -Ci) 



where := X'(s) 9 ^x(s)) ■ That implies 



(7.48) 



(7.49) 



(7.50) 



and 



- 2w'oo = ^ a^Zi^Uoo - ei) + + ej)) 



(7.51) 



i.e. 



X(u) = ^X'(u) ^ ai[2Z(u - d) + Ziuoo + e,) - ^(woo - e*)] 



(7.52) 



We also have: 



y{u) = -y\u)Y,^z{u 



d) + Z(Moo + d) - Z(U C 



(j>(s — u) 
X'{u) 



(7.53) 



In particular near u = e*, we find: 

1 X"'(e) 

and near tt = Moo, we write: 

#(u) (7.55) 

U - Moo 

— = a4{u - Uoo) (7.56) 
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and near u = — Woo, we write: 

y(u) ^— (7.57) 

^ = ^2a i (f)(u + u 00 ) (7.58) 

% 

note that: 

A-S^ A--~S^- (7 59) 

Note also that: 
7.0.2 variation of y(e;) 



-^(ei) ^a^fe - e,) - Ci^j + ^y"^) 
1 v „ \ ^ w (e,) 0(e,- S ) /l^(e,) 1 A"" 2 (e 



+ 



(7.61) 



2 v ,y A"'(e 4 ) *"( ei ) V6 <*"'(e*) 4 A"' 2 (e;) 
l <t>'{e i -s)X'"{e i ) i V'jej-s) 

2 A"' 2 (e,) 2 X"(e t ) 



i.e. 



^(lny(e,)) ^ 



+ a^fy'"(e t ) t lX IV (e z ) IX'" 2 ^) y"( ei )X'"(e t 



2 V y'ifr) 3 ^"(e<) 2 A"' 2 (e,) y'(e<) *"(e 

1 (j/jej - s)X'"( ei ) _ 1 (j>"{ej-s) 

2 A"' 2 ( ei )y(e,) 2*"(e i );y(e i ) 

(7.62) 

thus: 



a, /y"( ei ) , l^( ei ) lA"" 2 ! 



2 V y(ei) 3 A"'(e,) 2 A"' 2 (e, 
y'fe) * w (e.) 
y (eO *"( ei ) 
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1 \- <t>'{ej - s)X"\e l ) 1 ^ ( j ) "(e l - s) 



X" 2 (e, i )y'(e l ) 2^X"(e i )y(e i ) 
(7.63) 

Notice that: 

= Z ( u - e *) - 2Z ( U - M -) - (d 2 + l)Z{u + Moo ) (7.64) 
implies (expand to order 1 in u — ej): 



therefore 



v Oi ( y'"(e t ) X IV ( ei ) X"'\e 



- 2 V y(ei) X"( ei ) X"*(e t 

y»(e t ) X"'(e t ) 

y(e t ) X"{ ei ) 
I 0^ - s )X'"{e i ) _l sr J%e t -s) 



A"' 2 ^)^) 2 ^ X"( ei )y(e t ) 



(7.66) 



i.e. 



v Oi f y"( ei ) X IV ( ei ) X'"\ ei ) 



2 V^(e,) X"( ei ) X"\e,, 

Tjej) X"'(e t ) 

y(e t ) X"{ ei ) 
1 ^ <f>'(ej - s)X"'{e i ) 1 ^ 0"(e, - s) 



2^ X'^y^) 2 ^ X"( ei )y (eO 



(7.67) 

therefore: 

24 

where -K" may depend on V2 and e but not on Vi. Using: 



= K(V 2 , e) - In #'(0) 8 A 2 i (rf2+1) JJ y (e,) (7.68) 
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we can write it in a more symmetric form: 



24 v ' 11 6»(e; - u 00 )6'(e i + u c 



(7.70) 



where .fT depends neither on V\ nor on V2. K may still depend on e. 

The calculation of derivatives wrt e is exactly the same, with replacing a, by 
1/X"(ei)y'(ei), and we find that .ff does not depend on e. 

The limit where e — > is the genus zero case, and we have X = 0. 

8 Conclusion 

We have computed the free energy to order 1/N 2 , in the genus one case, and it should 
not be too difficult to do the calculation for all genus. It is conjectured that F^ 1 ' should 
be the log of the determinant of the Laplacian on the algebraic curve. 

Aknowledgements: I am thankfull to M. Bertola, V. Kazakov, I. Kostov for fruitful 
discussions. 

Appendix A Calculation 

We shall prove here that: 

Notice that the second of these two equalities follows from the first. 
The compatibility condition coming from eq. ([2.17)1 : 

d dF<® d dF<® d rir , N d -*r/ \ . . 

-W x {x) = -—Y(x) (A.2) 



dVi(x) dej dej dVx(x) dej de 3 

as well as eq. (j3.17j) and eq. f)3. 19j) imply that — is independent of Vi, and by 
the same argument after an integration by parts, it must be independent of V 2 too. 

In order to prove that T 3 -^ = 0, we can choose V\ and V2 such that the x- 
physical sheet contains all the A-cylces. That means that the function Yq(x) has g + 1 
cuts in the physical sheet, i.e. the large average density of eigenvalues of matrix Mi 
has a support made of g + 1 connected intervals: 

9+1 

'Y(x 4- iO) — Y(x — iO)) . sunn o = 

2ztt 



I 9+1 
p(x) = —(Y{x + iO)-Y(x-iO)) , Buppp=|J[a i ,6 < ] (A.3) 



=1 
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We choose the A contours as follows: 




and we can choose V\ and V2 such that the filling fraction of interval [a*, 6j] is: 



N 



p(x)dx 



1 

2nr 



yds = e 4 - 6i_i 



(A.4) 



a,. 



rij is the average number of eigenvalues of Mi in interval [a^, bi\. 

Now, following the method of for the one-matrix model, compute the energy 
cost of moving one single eigenvalue x from a to b, where a belongs to cut i and b to 
cut i + 1. Basicaly: rij — > rii — 1 and n.j+i — > rij+i + 1, i.e. e, — > — 1/iV. That energy 
cost can be written: 

1 dN 2 F f b 

-N—= N L w,dx (A ' 5) 

where the effective mean field potential V e g experienced by one eigenvalue of M\ in 
the presence of the potentials and the interaction with all other eigenvalues in their 
equilibrium position was computed in , and it reduces to: 



(A.6) 



BF 1 

_ = —eh Y{x)dx = r< 



dei 1m j Bi 



Appendix B Theta functions in genus one 

Consider a complex number r such that Imr > 0, called the modulus. We define the 
theta-function by: 



Q(u) = J2e inrn2 sin (2™) 

71=0 

We have: 

6{u + 1) = 9{u) , 6{u + r) = -e{u)e- i7T ^ u+ ^ 
We also introduce: 



we have: 

Z{u + 1) = Z(u) 



Z{u) = ^-\n9{u) 



Z(u + t) = Z(u)-2i7r 



(B.l) 

0(- u ) = -6(u) (B.2) 
(B.3) 

Z(-u) = -Z(u) (B.4) 
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Z(u) has a single pole at u — 0, with residue 1: 



~ - + Ci« + 0(m 3 ) (B.5) 
And we introduce the Weierstrass function: 

= (B.6) 

we have: 

0(w + l) = 0( m + t) = <f>(-u) = <j)(u) (B.7) 
is thus an elliptical function (doubly periodic). It has a double pole at u — 0: 

1 



0H~^-G + O(« 2 ) (B.8) 
Consider an elliptical function / such that: 

/(« + 1) = /(« + r) = /(«) (B.9) 

then: 

• if / is entire, then / is a constant. 

• if / is meromorphic (its singularities are poles), it must have at least two poles 
(possibly a double pole). 

• if / has n poles ei,...,e n with multiplicities m 1? ...,m n , in the fundamental 
parallelogram of side (1,t), then there exist m = X]fc=i m & complex numbers 
A, /i, f 2 , . . . , f m -i such that: 

f( U ) = ^y^" 1 ?' (b.io) 

where we have defined f m such that: 

m n 

"}2fk = ^2rn k e k (B.H) 

k=l k=l 

• An alternative representation of / is the following: there exist there exist m + 1 
numbers A , A hl , A hmi , A nA , A U)mn such that: 

n m k 

f(u) = A + J2Y, - e fc ) (B.12) 

fc=i z=i 

with the condition that: 

n 

J]A M = (B.13) 

fc=i 
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Appendix C Theta functions arbitrary genus 



C.l Genus g > 1 

e( u ) = J2e inntrn e 2inntu (C.l) 

n 

we have: 

0( u + n) = 6{u) = 6{-u) (C.2) 

and 

Q{u + rn) = 6{u) e ^2«W™] ^ ^ 
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